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We consider spin Hall effect in a system of massless Dirac fermions in a graphene lattice. Two 
types of spin-orbit interaction, pertinent to the graphene lattice, are taken into account - the intrinsic 
and Rashba terms. Assuming perfect crystal lattice, we calculate the topological contribution to 
spin Hall conductivity. When both interactions are present, their interplay is shown to lead to some 
peculiarities in the dependence of spin Hall conductivity on the Fermi level. 
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I. INTRODUCTION 

Long time ago Dyakonov and Perel predicted, that in 
some systems transverse spin current and spin accumu- 
lation may be induced by electric current 1 ^. This effect, 
known as spin Hall effect (SHE), is currently extensively 
studied both experimentally and theoretically^i^&ii^ 
(see also review paper o 9 ' 10 ). The effect may appear in 
semiconductors as well as in metals, and originates from 
spin-orbit interaction. Generally, such an interaction 
may be of intrinsic or extrinsic origin. The extrinsic SHE 
is associated with scattering mechanisms, like skew scat- 
tering and side jump in the presence of impurities. In 
turn, the intrinsic mechanism of SHE is a consequence of 
an unusual trajectory of the charge carriers in the mo- 
mentum space, which may be described by the Berry 
phase formalis m 11 ! 12 . This contribution will be referred 
to as the topological one. 

In this paper we consider the topological contribution 
to spin Hall conductivity in a system of Dirac fermions 
with spin-orbit coupling. Generally, the form of spin- 
orbit interaction depends on the symmetry and structure 
of the system. The Dirac model turned out to be useful 
not only in the relativistic field theory, but also in con- 
densed matter physics to describe some features of elec- 
tronic spectrum (at least in a certain energy range). One 
of such systems is two-dimensional graphene, and in this 
paper we consider Dirac fermions with the spin-orbit in- 
teraction taken in the form appropriate for graphene and 
including both intrinsic and Rashba terms. We believe, 
that the results derived here will shed some light on the 
spin Hall effect in graphene. 

Graphene is a two-dimensional honeycomb lattice of 
carbon atoms, with two noncquivalcnt sublattices. The 
low-energy electron states near the K and K' points at 
opposite corners of the Brillouinc zone can be approx- 
imated by the conical energy spectrum. As a result, 
charge carriers are described by the Dirac equation 13 i 14 . 
The Fermi surface in a neutral graphene consists of the 



nonequivalent points K and K', at which the valence and 
conduction bands touch each other. However, when the 
intrinsic spin-orbit interaction is included, an energy gap 
opens at these points. Unfortunately, it is now believed 
that the intrinsic spin-orbit coupling in graphene is rather 
weak so the gap is also small. Since the graphene layer 
is usually on a substrate, one also can expect spin-orbit 
interaction of Rashba typei^. Moreover, the correspond- 
ing coupling parameter can be tuned externally by a gate 
voltage. Indeed, a large Rashba spin-orbit interaction has 
been reported in a recent experiment — . 

Kane and Meleil have shown that the intrinsic spin- 
orbit interaction opens an energy gap at the Dirac points, 
and also have predicted a quantized value of the spin 
Hall conductivity when the Fermi level is in the gap. 
The quantized spin Hall conductivity at the Dirac points 
was also confirmed by later analytical and numerical cal- 
culation o 18 i 19 ' 20 . On the other hand, the presence of 
Rashba spin-orbit interaction reduces the gap, and when 
the Rashba interaction is stronger than the intrinsic one, 
the gap becomes closed. Since the magnitude of intrin- 
sic spin-orbit interaction seems to be significantly smaller 
than that assumed originally—, it is possible to reach the 
limit opposite to that considered by Kane and Mele, i.e. 
the limit where the Rashba coupling dominates while in- 
trinsic spin-orbit interaction is negligible. To our knowl- 
edge, spin Hall effect in this limit has not been considered 
analytically so far. Therefore, in this paper we present 
analytical and numerical results obtained within the lin- 
ear response theory and Green functions technique, as- 
suming both intrinsic and Rashba spin-orbit interaction. 
We focus on the topological contribution to the effect 
assuming perfect crystal lattice and ignoring impurities 
and defects, which however may influence the magnitude 
of spin Hall effect in real system a 18 ' 19 ! 20 ' 22 . 

The paper is organized as follows. In section 2 we de- 
scribe the model and present a general formula for the 
spin Hall conductivity. The case of pure intrinsic spin- 
orbit interaction is presented in section 3, while the case 
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of Rasba interaction is described and discussed in section 
4. The general case, where both interactions are present, 
is described in section 5, while summary and final con- 
clusions are in section 6. 



II. MODEL AND GENERAL FORMULA FOR 
SPIN HALL CONDUCTIVITY 

Including both intrinsic and Rashba spin-orbit cou- 
pling, the effective mass Hamiltonian 1 ? of graphene can 
be written in the form 



with / being the 2x2 unit matrix. 

In the linear response theory, the dc spin Hall conduc- 
tivity is given by the formula^, 



lim — Tr 



de d 2 k 
2^ (2tt) 2 



[vx,8 z ] + G k (e+w)u y G k (e), 

where Gk(e) is the causal Green function corresponding 
to the Hamiltonian (1). This formula will be used in 
the following to calculate spin Hall conductivity in some 
specific cases as well as in a general situation. 



H = H n 



H. 
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The first term, Hq, describes the low-energy electronic 
states around the Dirac points K and K' in the Brillouin 
zone, and has the form 



H 





v(±k x + ik v ) 



v(±k x - iky) 




(2) 



where the upper and lower signs correspond to the points 
K and K', respectively, and v is a parameter describing 
the conical energy spectrum, v = Hvf, with vf denoting 
the electron Fermi velocity. The second term in Eq.(l) 
describes the intrinsic spin-orbit interaction in graphene, 



H SO = 



±A so <t z 
tAsoo-* 



(3) 



with Ago being the relevant parameter (2 Ago is the gap 
created by the intrinsic spin-orbit coupling in the Dirac 
points). Finally, the last term in Eq.(l) stands for the 
Rashba spin-orbit term, 



Hp. 







A/e(±CT a -I- ia x ) 




(4) 



where \r is the corresponding coupling parameter. As 
in Eq.(2), the upper and lower signs in (3) and (4) cor- 
respond to the two inequivalent points K and K' of the 
Brillouin zone, respectively. 

To obtain spin Hall conductivity we introduce first the 
spin current density operator, 



f 3 = \ [v, Sj] + 



(5) 



III. THE CASE OF X R = AND A SO + 

We will consider first the special case, when the Rashba 
coupling vanishes, while the intrinsic spin-orbit interac- 
tion is nonzero, A#, = and Ago 7^ 0. Such a situa- 
tion has been already considered analytically in the clean 
limitii, and also studied numerically in the presence of 
impurities^. Our results are consistent with those ob- 
tained in the above cited works. 

When Xji = 0, Eq.(8) for the spin Hall conductivity 
takes the following form: 
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where /i is the chemical potential. 

In the zero temperature limit, the spin Hall conductiv- 
ity will be presented in the form 



xy 



(10) 



where <J X y'° is the contribution from the fully occupied va- 
lence band, and Sa xy is associated either with the empty 
part of the valence band (upper sign) or with occupied 
part of the conduction band (lower sign) . 



The contribution 5a x z y takes the form 
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where [A, B}+ = AB + BA denotes the anticommutator 
of any operators A and B, Vi = (\/Ti)(dH / dki) is the 
velocity operator (i — x,y) and Sj — (%/2)aj is the j-th 
component (J = x, y, z) of the spin operator. Taking into 
account the exact form of Hamiltonian H one finds 



= v I 

h \i 



and 
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where the upper sign refers to the situation when the 
Fermi level is in the valence band, while the lower one 
when the Fermi level is in the conduction band. In turn, 
the contribution from the fully occupied valence band 
can be calculated as 
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FIG. 1: (color online) Spin Hall conductivity in the absence 
of Rashba interaction and for indicated values of the intrinsic 
spin-orbit parameter Aso- The inset shows the energy spec- 
trum in the vicinity of the Dirac point for Aso = O.OlmeV. 
The parameter v is defined as v = % Vf, where vf is the Fermi 
velocity, v F = 0.833 x 10 6 m/s. 



Thus, from Eqs ([TTJ and |TJ} follows that the whole 
spin Hall conductivity can be written as 



*so 



for \/j,\ > Aso, and 
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(14) 



for \/i\ < Aso- The above result coincides with the ones 
obtained by Kane and Meleii and Sinitsyn et ali2. The 
other two component of the spin Hall conductivity, i.e. 
a%y and a X y vanish, as one could expect. 

The corresponding numerical results are shown in Fig.l 
for three different values of the parameter Aso- Addi- 
tional factor of 2 has been taken into account in order 
to include the two Dirac points K and K'. The spin Hall 
conductivity is shown there as a function of the chemical 
potential measured from the middle of the gap. When 
the Fermi level is in the energy gap, the Hall conductiv- 
ity is constant and quantized. When, in turn, the Fermi 
level is either in the valence or conduction bands, the ab- 
solute value of conductivity decreases, and disappears for 
\fi\ — » oo. Note, the conductivity is symmetric with re- 
spect to the middle of the gap. The inset in Fig.l shows 
the energy spectrum in the vicinity of the Dirac point, 
calculated for the parameter vp obtained by Gmitra et 

„124 



IV. THE CASE OF A SO = AND \ R / 

Now, we consider the opposite situation, i.e. when the 
intrinsic spin-orbit coupling is negligible, Aso — 0, while 



the Rashba parameter is nonzero, X R ^ 0. Equation ((SJ) 
leads then to the following formula for the spin Hall con- 
ductivity: 

de rf 2 k 8iev 2 X%(e + fi)[v 2 (kl-kl) + {e + fi) 2 

^ ( 2 ^) 2 ULi [* ~ E n (k) + y. + iS sign(e)] 

1 

x 4 ■ 

n m =i[ £ - £»( k ) + M + ^sign(e)] 
1 — 4) describe the electron energy 



where Ei(k) (i 
spectrum, 



Ei (k) = A/ 

E 2 (k) = X R -(X R - 
E 3 (k) = -X R + (X 2 R - 
E 4 (k) = -X R - (A| 



Vfc 2 ) 1/2 , 

v 2 k 2 )^ 2 , 
v 2 k 2 f'\ 

-v 2 k 2 Y' 2 . 



(16) 
(17) 
(18) 
(19) 



The states E\ (k) and E 3 (k) correspond to the conduction 
bands while ^(k) and ^(k) to the valence bands. We 
consider first the case when > 2 X R . 



A. The case of \fj,\ > 2 Xr 

When the Fermi level is in the two valence bands, /1 < 
—2X R , then upon integrating Eq. (|15p over e one arrives 
at the following formula 



167rAi?, 



d k 2X » k 



(v 2 k 2 + A|) 3 / 2 



if(e 2 )-f(e 4 )}, (20) 



where /(e) is the Fermi-Dirac distribution (assumed here 
for T = 0). Taking now into account the notation intro- 
duced in Eq.(10), one finds 



0. 



and 
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(21) 



(22) 



Att 2( M 2 - A|) ■ 

As before, a^y is the contribution from both fully oc- 
cupied valence bands, while Sa^ z y takes into account the 
empty part of the valence bands. Now, the contributions 
from the two fully occupied valence bands cancel each 
other, so a X y° vanishes exactly. 

When the Fermi level is in both conduction bands, /1 > 
2 X R , the spin Hall conductivity can be calculated in a 
similar way, and the formula for Sa xy takes the same form 
as for [i < 2 X R , i.e. Eq.(22). Thus, the total spin Hall 
conductivity for > 2 X R can be written as 



(23) 



xy 2{fi 2 - X 2 R ) 4tt 
for the Fermi level in conduction bands, fi > 2 X R , and 
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2{y? - X 2 R ) 4tt 
for the Fermi level in the valence bands, [i < 2 X R . 
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FIG. 2: (color online) Spin Hall conductivity in the absence of 
intrinsic spin-orbit interaction and for indicated values of the 
Rashba coupling parameter Xr. The other parameters are as 
in Fig.l. The inset shows the energy spectrum in the vicinity 
of the energy gap for X R = O.OlmeV. 



B. The case of < 2 X R 

Now, we assume that |^| < 2 Xr. The only difference 
is that now the Fermi level is either in one valence band 
or in one conduction band. Spin Hall conductivity can 
be calculated in a similar way as before and one finds 



fi(^ + 2X R ) e 



'xy 



4X R (n + X R ) 4tt 
for the Fermi level in the conduction band, and 



(25) 



'xy 



Mm - 2 \r) 

4Afl(/i - X R ) 4tt 



for the Fermi level in the valence band. 

Figure 2 presents numerical results for the spin Hall 
conductivity as a function of the chemical potential. 
As before, we included a factor of 2 in order to take 
into account contribution from the second Dirac point. 
We note, that now the spin Hall conductivity tends to 
a xy — ^ e/47r in the limit of /i — » — oo, while for [i — > oo 
it tends to the a x z y = e/4ir. This behavior is different 
from that for intrinsic spin-orbit interaction, where the 
contribution from conduction bands cancelled the corre- 
sponding part from the valence bands. It is also worth 
to note that now the spin Hall conductivity is antisym- 
metric with respect to change of the Fermi level sign. 

We have also checked the other components of the spin 
Hall conductivity. As before, these components vanish 
exactly, a s x x y = o x 



'xy — 0. 



V. THE CASE WITH X R / AND A so + 

When both intrinsic and Rashba spin-orbit interac- 
tions are present in the system, their interplay leads to 



interesting and peculiar behavior of the spin Hall con- 
ductivity. The analytical formulas, however, are much 
to complex to be presented here, so we will show mainly 
results of numerical calculations. 
Writing the Green function as 



G k (e) 



5k(e) 



II,Ui[ e - ^n( k ) + A* + itfsign(e)] 



(27) 



where E n (k) (n = 1 — 4) are the dispersion relations for 
the conduction and valence bands, 



£i(k) = X R + ^(A so -\ R ) 2 + v 2 k^ 
E 2 (k) = X R - V(Aso-A«) 2 + « 2 fc2, 
E 3 (k) 
E 4 (k) 



X R + ^{A so + \ R ) 2 +v^k\ 
Xr 



V(A50 + A fl ) 2 + « 2 fc 2 



(28) 
(29) 
(30) 
(31) 



and taking into account the explicit form of v x and v y 
one may write the relevant trace in Eq.(8) as 



Tr 



a z 

G Z 



9k ( £ i 



-/ 
/ 



9 k (e 2 )\ =V(e 1 ,e 2 ^2) 



where V{e\, e 2 ) is a certain function of E\ = e + fj, ■ 
and e 2 = e + /i. Expanding V(e\, e 2 ) as 



V(ei, e 2 ) = V(ei, E 2 )| El=E2=e+M 

dV(ei, e 2 ) 



+0J 



and taking into account that V(si, s 2 ) 
one finds 



£l=£2=£ + V 



..(33) 
0, 
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+16A 2 ? ( £ + /i) 3 -4A so (e + M) 4 }- (34) 



These formulas can be then used to calculate spin Hall 
conductivity, either analytically by integrating over e and 
k or numerically. Since analytical formula are gener- 
ally rather cumbersome, we performed numerical calcu- 
lations, while analytical formula will be presented only 
for some special cases (see below). 



A. X R > A so 

Consider first the case when is significantly larger 
than Aso- The corresponding spin Hall conductivity 
is shown in Fig. 3 for X R = 0.1 meV and three differ- 
ent values of Aso (smaller than Xr). When Aso is 
much smaller than A^, then the conductivity is deter- 
mined practically only by the Rashba coupling. With in- 
creasing Aso ! the interplay of both interactions leads to 
anomalous behavior of the spin Hall conductivity. More 
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FIG. 3: (color online) Spin Hall conductivity for A_r = O.lmeV 
and indicated values of Aso, A so < Xr. The inset shows 
the energy spectrum in the vicinity of the Dirac point for 
Aso = O.OlmeV. The parameter v is taken as in Fig.l. 
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FIG. 4: (color online) Spin Hall conductivity for Xr = 
O.OlmeV and indicated values of Aso, Aso > Xr. The inset 
shows the energy spectrum in the vicinity of the Dirac point 
for Aso = 0.05meV. The parameter v is taken as in Fig.l. 



specifically, the conductivity becomes diverging when the 
Fermi level approaches the point, at which the top of 
the upper valence band touches the bottom of the lower 
conduction band (see the inset in Fig. 3). This appears 
when fi — Aso- Asymptotic behavior of the conduc- 
tance near the point /j, = Ago is described by the 
term (A so /2A_r) ln[-A so - Xr + \/ (Xr + m) 2 ] when 
/i = Aso is approached from the right (/i > Aso) side, 
and (A so /2A fl ) ln[A so - X R + \/{X R ~ n) 2 } wn en it is 
approached from the left (/i < A S o) side. 



comparable. Some simple analytical results, however, 
can be obtained for Xr = Aso- For this particular case, 
the bottom edges of two conduction bands coincide with 
top edges of one of the valence band, while the top edge 
of the second valence band is much below (see the inset 
in Fig. 5). 

The relevant formula for the spin Hall conductivity 
depends then on the Fermi level as follows: 



1. fi< -3Aso 



B. 



Xr < Aso 



Let us now consider the opposite situation, when Xr is 
smaller than Aso- The corresponding spin Hall conduc- 
tivity is shown in Fig. 4 for Xr = 0.01 meV and three dif- 
ferent values of Aso (larger than Xr). General shape of 
the curve showing spin Hall conductivity as a function of 
the chemical potential is similar to that for Xr = 0. How- 
ever, the interplay of intrinsic spin-orbit interaction and 
Rashba coupling leads to an interesting feature. More 
specifically, there is now no divergence, but a weak kink 
in the conductance appears on the negative chemical po- 
tential side. When Aso decreases and approaches A^, 
the kink becomes more pronounced. The kink is asso- 
ciated with splitting of the valence band by the Rashba 
interaction. The upper valence band edges are now at 
-A so - 2Aj? and -A S o + 2A i? . 



Ac 



Variation of the spin Hall conductivity with the Fermi 
level becomes more complex when both Aso arL d Xr are 



For \i < — 3Aso, the spin Hall conductivity is given 
by the formula 



s z _ _ 

°xy 



±so 



■in 



(35) 



This formula covers the energy range up to the top edge 
of the lower valence band. 



2. 



-3Aso < fJ, < Aso 



When —3 Aso < M < Aso, the corresponding formula 
for the conductivity takes the form 



7 S = 
'xy 
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so- V 
4A 



so 
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(36) 



This formula, in turn, describes spin Hall conductivity 
when the chemical potential is between the top edge of 
the lower valence band and bottom edges of the conduc- 
tion bands. Note, that this formula leads to diverging 
spin Hall conductivity when the Fermi level tends from 
left (lower values) to Aso- 
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3. n> A so 

Finally, for fi > Aso the spin Hall conductivity is 
equal to 

a - = VTa^^ (37) 

This formula gives a finite spin Hall conductivity in the 
whole range of its applicability, also at the point /i = 
Aso- 




\jl (meV) 

FIG. 5: (color online) Spin Hall conductivity for indicated 
values of \r = Aso- The inset shows the energy spectrum 
in the vicinity of the Dirac point for A_r = Aso = O.OlmeV. 
The parameter v is taken as in Fig.l. The vertical dotted 
lines indicate the position where the anomaly appears. 

The spin Hall conductivity for Xr — Aso is shown in 
Fig. 5 for indicated values of Xr = Aso- The anomaly 
at jj, = Xr = Aso is now clearly visible. When the 



chemical potential y, tends to /i = Xr = Aso from the 
right (n > Aso) side, the conductivity is finite, while 
when it tends to fi = Xr = Aso from the left (/x; Aso) 
side, the conductivity becomes diverging. The vertical 
dotted lines in Fig. 5 indicate only the position where 
the anomaly appears. 

VI. SUMMARY AND DISCUSSION 

Assuming intrinsic and Rashba spin-orbit interaction 
we have calculated topological contribution to the spin 
Hall conductivity. In the limit of vanishing Rashba term 
we arrived at the results which are in agreement with 
those available in the relevant literature. When, in turn, 
the Rashba coupling dominates and intrinsic spin orbit- 
coupling vanishes, we have found asymmetric behavior 
of the spin Hall conductivity with respect to the sign 
reversal of the chemical potential. Such a change in the 
chemical potential can be achieved with a gate voltage, 
for instance. 

When both intrinsic and Rashba spin-orbit interac- 
tions are present, their interplay leads to some pecu- 
liar and anomalous behavior of the spin hall conductiv- 
ity with the Fermi level. In particular, for some range 
of spin-orbit parameters, the spin Hall conductivity was 
found to diverge when the Fermi level approaches the 
limit fi = Aso- 
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